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A. Voigt, J. Phys. A: Math. Gen. 27, 1139 (1994) [1]. The model is defined by 
H = >/i X^lili ^oSj- + J2HR{si} ; Ji, J2 > , i = 1,...,N. In extension to Ref. [1] 
we consider a more general s,} describing the properties of the spins surrounding 
the central spin sq. The Heisenberg star may be considered as an essential structure 
element of a lattice with frustration (namely a spin embedded in a magnetic matrix 
Hr) or, alternatively, as a magnetic system Hr with a perturbation by an extra spin. 
We present some general features of the eigenvalues, the eigenfunctions as well as the 
spin correlation (sqSj) of the model. For Hr being a linear chain, a square lattice 
or a Lieb-Mattis type system we present the ground state properties of the model in 
dependence on the frustration parameter a = J2/ Ji- Furthermore the thermodynamic 
properties are calculated for Hr being a Lieb-Mattis antiferromagnet. 
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1. Introduction 

The properties of interacting quantum spins have been attracted large attention over a 
long period. Only some model Hamiltonians can be solved exactly. Important examples 
are (i) models in one dimension solvable by Bethe-Ansatz [2, 3, 4], (ii) valence-bond 
models [5, 6, 7], (iii) a one-dimensional model with long-range inverse-square exchange 
[8, 9, 10, 11], and (iv) models with long-range interaction of constant strength (Lieb- 
Mattis type models [12] which have been used to discuss spontaneous symmetry breaking 
in quantum spin systems recently [13, 14, 15]) . 

In [1] (further referred as STAR I) we introduced the frustrated spin | Heisenberg 
star with a Hamiltonian Hi = Ji Ylf=i ^o^i + J2 Ylf=i ^i^i+i {Ji^ J2 > 0) representing a 
central site with N nearest neighbours which can be unconnected ( J2 = 0) or connected 
{J2 7^ 0). The star can be considered either as an essential structure element of a lattice 
or as a linear chain with a perturbing extra spin. 

For the above defined Hamiltonian Hj we presented in STAR I general relations for 
the eigenvalues, the eigenstates and the spin correlation function in the ground state as 
well as numerical results for N = 4, 6,..., 22. Analyzing the analytical and numerical 
data we discussed the ground state phase diagram, in particular, the ground state spin 
correlations versus a = J2/ Ji- We found that for J2/ Ji < ctcrit the ground state of the 
system is the state with strongest antiferromagnetic correlation (soSj) = — between 
the central spin Sq and a neighbouring spin and with ferromagnetic correlations 
(Sj-Sj) = I within the ring. If J2/J1 exceeds acrit it follows a series of transitions to 
states with successively weaker correlations (soSj) ending with (soSj) = for dominating 
J2. For acrit we found exactly | , independent of the size of the system. For larger N this 
weakening of the antiferromagnetic correlation of the central spin takes place very rapidly 
when changing J2/ Ji in a small region above acrit- The extrapolation N ^ 00 yields 
for J2/J1 < acrit the correlator (soSj) = — | which can be considered as an upper limit 
for the ground state correlation (s;Sm)o of antiferromagnetically interacting spins / and 
m in a spin | Heisenberg antiferromagnet without competition between the interactions. 
We argued that any ground state spin correlation (sjSj) of antiferromagnetically coupled 
spins Si and Sj larger than — | is an effect of competing interactions. In the limit of large 
J2I Ji the correlation {siSj)o within the ring becomes antiferromagnetic and the ring state 
goes over to the Bethe singlet for J2 ~ JiN{-^ -\- 0{j^)). 

In this paper we generalize the model by considering a central spin with different 
embedding media. In particular, we compare models where a central spin is embedded 
in an antiferromagnetic Heisenberg matrix of linear chain, square lattice and Lieb- 
Mattis type. Furthermore, we discuss in more detail the properties of the model in 
the thermodynamic limit N ^ 00. Additional to the ground state properties we present 
the full thermodynamics in the case of the Lieb-Mattis star. 
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2. Model and general relations 



We consider a spin | Heisenberg system 



J I ^ 

H = —J2^oSi + J2Hr{s,} ; Ji,J2>0 

8 = 1 



(1) 



H 



R 




Here HR{si} represents itself a spin | rotationally invariant antiferromagnet describing 
the medium which surrounds the central spin. In (1) the interaction with the central spin 
was scaled by N. A similar scaling will also be used for Hr (see below). In contrast to 
the paper STAR I we do not specify the Hamiltonian Hr at this point. 

Several general relations found for the LC star in paper STAR I are valid for the 
more general Hr^ too. Let us define the total spin S = YliLo of the system and the 
total spin of the embedding medium S^^ = S — Sq = Y^f-i s^. The following universal 
features of the model are important: 

(i) The integrals of motion of the systems are i7, Hr^ S^, S*^, S|j with the respective 
quantum numbers Er^ s, m, r. 

(ii) The eigenvalues of H are given by 

E = .hER + ^r for 3 = r + ^ ; r = 0, 1, 2, | (2) 

E = J,ER-^{r + l) for 3 = r - ^ ; r = 1, 2, |. (3) 

Clearly, for fixed r and fixed Er the states with r = s + 1 have lower energy. Let us look 
for the lowest energy in every subspace of fixed quantum number r. For the energy Er of 
the embedding antiferromagnet Hr we have the Lieb-Mattis level ordering [12, 16] for the 
lowest eigenvalue ER^r) in the considered subspace. Hence, from the general relation (3) 
we indicate a competition between both energy terms leading to interesting frustration 
effects. In particular, at zero temperature the frustration parameter a = J2/ Ji determ- 
ines the quantum number r in the ground state and may be used to tune the ground state 
properties. 
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(iii) The basic structure of the eigenstates of (1) reads as follows 

\^E,s,m,r) = a I t )l$L.,™-l/2) + H i ) I $L,.,™ + l/2) , (4) 

where | '\ ) and | J, ) are the eigenfunctions of the z-component of the central spin So,^ 
and the \^^^r m±il2) eigenfunctions of Hr^ S^, Sr^z with the corresponding 

eigenvalues Er, r{r + 1), m ± 1/2. The eigenfunctions |$i^,r,m+i/2) ^^^^ |$i^,r,m-i/2) 
of i^i^ are related to each other by 

l$L.,™-i/2) = [r{r + 1) - (m + l/2)(m - 1 12)]-''' SrM^^,^^^,,,) (5) 
The coefficients a and h do not depend on Ji, J2 and are given by 



r+m+l r — m + ^ ^ 1 



2r + l V 2r + 1 2 

/r — m + l /r + m + l ^ 1 



for s = r H — (6) 



, . ,, for s = r . (7) 

V 2r + 1 ' V 2r + 1 2 ^ ' 

(iv) For the spin correlation of the central spin with a spin of embedding medium we 
have 

r 1 

{^E,s,m,r\^0Si\^E,s,m,r) = _ for S = T + - (8) 

-(r + 1) 1 

{^E,s,mA^O^^\^E,s,m,r) = _ for S = T - - . (9) 

The spin correlation within the embedding medium, of course, depends on the detailed 
structure of Hr. 

The more specific properties of the model are determined by properties of the 
embedding medium entering the general relations (2) - (9) via Er, \^ERrm±i/2)- 
Whenever the solution for the Hamiltonian of the surrounding medium Hr is known 
(i.e. Er and I^Ert ■m±i/2) available) the total system is solvable. 

We consider as solvable systems the linear chain (LC) and the Lieb-Mattis 
(LM) antiferromagnet. Both model represent the extreme limits of the Heisenberg 
antiferromagnet with maximal (LC) and minimal (LM) quantum fiuctuations. 
Additionally to the solvable limits we present approximative results for the square-lattice 
(SL) antiferromagnet. In order to compare the different cases we scale the interactions 
by the number of sites N so that the energy becomes intensive and that the energy of 
the fully polarized ferromagnetic state is identical E{^ = | for all embedding media 
considered. This scaling is consistent with the 1/A^-scaling of the interaction of the 
central spin (see equation (1)). We write for the Hamiltonians of the media 

^^"^ =tfE«»-«hi (10) 

8 = 1 
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rSL 

R 



2N 



1 




(11) 
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iV2 



N 



rLM 

R 



(12) 



In what follows we fix Ji = 1 and consider a = J2/ Ji as the parameter of the model. 



3. Finite systems 

For small N the star can serve as an elementary cluster of a lattice with frustration. 
We focus our consideration in this section on the example of N=8. For that N the LM 
system is by chance identical with the SL. For comparison we add to the discussion of the 
LC, SL and LM star the elementary cube of the body-centered cubic lattice (BCC). The 
calculation of the spectra and the wave function can be simply done by numerically exact 
diagonalization. The results for the relevant spin correlation of the central spin Sq with 
a neighbouring spin s^ are presented in figure 1. The steps in the correlation functions 
indicate the transitions between states with different quantum number r. Of course, the 
frustrating J2 interaction weakens the antiferromagnetic correlation. But it is interesting 
that independent of the embedding medium an essential diminishing of the strength of 
the correlation takes place in a small region above acrit = {J2I Ji)crit = \- This region is 
in particular small if the surrounding medium has a low number of nearest neighbours 
(LC). Otherwise, the complete suppression of the antiferromagnetic correlation by the 
frustrating J2 takes place for fairly large a (precisely at for LM and at slightly lower 
values for LC and BCC). 

For small N one can calculate the full thermodynamics. As an example we present the 
temperature dependence of (soSj) and (sjSj) for the LC star in figure 2. (The behaviour 
is qualitatively the same for the other systems.) For a < | the behaviour is the standard 
one for both (soSj) and (sjSj), i.e. the strength of the correlation is diminished by thermal 
fiuctuations. Of course, there is no sharp transition for the finite system. For a > \ the 
frustration is more important and we find for low temperatures a qualitatively different 
behaviour for (soSj) and (sjSj). While the thermal fluctuations diminish the correlations 
within the medium at the same time they increase the strength of antiferromagnetic 
correlations of the central spin. This order from disorder phenomenon was observed in 
several frustrated systems [17, 18, 19, 20] and is connected with a competition between 
different energy scales. In particular, slightly above a = 2 the ground state correlation 
(soSi)o is zero (cf. flgure 1), but the fluctuations cause an antiferromagnetic alignment 
of the spins at flnite temperatures. 
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4. The thermodynamic limit 
4.1. Ground state 

Now we turn over to the thermodynamic hmit. Then the star represents an 
antiferromagnet frustrated by an interaction with an extra spin Sq. This situation 
is somewhat similar to the slightly doped high Tc cuprate superconductors where the 
holes at the oxygen sites create antiferromagnetically coupled extra-spins frustrating the 
antiferromagnetic copper matrix [21, 22, 23, 24]. 

We define the normalized quantum number of the medium 

2 . . 

x = — r ; 0<a;<l. (13) 

The ground state energy is obtained by equation (3) by selecting the lowest eigenvalue 
Er for a given quantum number r and finding that r which minimizes the energy for 
given Ji and J2. Er of the LM antiferromagnet is given by the analytic expression 

E]f="^- \. (14) 

For the LC no simple analytic expression is available, the £'|j'^(a;)-function was calculated 
by the numerical solution of the Bethe-Ansatz equations [3, 4]. However, near maximal 
polarization, x — > 1, we can extract from the Bethe-Ansatz an analytic relation for the 
energy Er 

Ejfix) = i - (1 - x) + J^(l - xf + 0((1 - xf). (15) 

For the square lattice only approximative results are known. We calculated E^{x) for 
lattices of = 16, 18, 20 and 24 sites with periodical boundary conditions. The energy 
scales with N~^^'^ [25, 26, 27]. By interpolation between the discrete points calculated 
for A^ = 16, 18, 20, 24 sites and extrapolation to A^ — > 00 we obtained the numerical data 
for E'^^ in the thermodynamic limit. The error of these SL data could be estimated 
by comparing our result for x = (antiferromagnetic singlet ground state) with best 
available results of various methods (variational Monte Carlo, world-line Monte Carlo, 
spin-wave theories, see the review [28]). We found an error of less then 3%. 

A general result (independent of the model and the size of the system) is the stability 
of the fully polarized [x = 1) ferromagnetic state of the embedding medium till precisely 
a = J. At this point a second order transition to a canted magnetic structure occurs. 
Only at infinite J2/J1 the antiferromagnetic singlet state is reached. In more detail 
we want to discuss the ground state energy, the spin-spin correlation and the order 
parameters. 

Energy - In figure 3 we present the ground state energy in dependence on a = J2/ Ji- 
We find a maximum slightly above a = |. For the LM star the maximum is at 
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a = I/a/S ~ 0.35. For the LC and the SL star the position of the maxima is obtained 
from the numerical data at 0.26 (LC) and 0.30 (SL). The maximum in E[x) indicates the 
region of strongest frustration and coincides according to the Hellmann-Feynman [29] 
theorem with the point where the nearest-neighbour correlation function (siSj)o vanishes 
(maximal spin canting). 

Spin-spin correlation - In figure 4 the spin-spin correlation functions (soSi)o and 
(siSj)o («, J 7^ 0, J = z -|- 1 for LC, j = i-\-xoTj = i-\-y for SL and j ^ i ^ B for 
LM) are presented. For the LM system these correlation functions are given by explicit 
formulae 

(sos^)o = - — — ; a > Q,c„t (16) 
ib a 

(s.-Sj)o = iTZ-^ - T ; a > a^„t ■ (17) 
32a^ 4 

For the LC and the SL stars they are a result of numerical calculations. It is evident 
that the physical properties are changed drastically in a small parameter region slightly 
above the critical acrit = \- For the LC system it is even suggested by the numerical 
data that the correlation changes at this point with an infinite slope. This can be checked 
analytically using the equation (15). We obtain with /3 = a — 

(sos.)o = - 1 + - + OifJ'^') ; a > (18) 

4 TT TT 

(s.s,)o= 1 - + + 0(/3^/') ; . (19) 

4 TT OTT 

In other words, for the infinite LC star the transition at a = | is extremely sharp but 
still of second order. The LM correlation has for /3 — > a finite slope and the correlation 
behaves according to (soSi)o = —\ + fS -\- 0{f3^) and (siSj)o = | — 4/3 -|- 0{f3^). The 
behaviour of the SL star is just intermediate between the two limiting cases LC and LM. 

Order parameters - On the base of the spin correlation we can define the order 
parameters 



1 JV 1 ^ 

i=l ^J = l 



for ferromagnetic and 

N 1 N 



Ms = {[J^J2W) = J^J2 r^r,{s,s,) ; r, = +1 , z G A ; = -1 , i e B [21] 

i=l ^J = l 



for antiferromagnetic long-range order within the medium. We choose here the definitions 
via the long-range part of the correlation functions and do not introduce symmetry 
breaking field (cf. [13, 14, 15]). The ferromagnetic order parameter can be expressed 
directly by the normalized quantum number x defined in (13), = jx^ and can be 
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calculated numerically for LC and SL, but analytically for LM. The antiferromagnetic 
order parameter can be evaluated for the LM system only. We have 

MIlM = T - ; « > (^crrt ■ (23) 

For the ferromagnetic order parameter M^q of the LC system an analytic expression is 
available near the second order transition, i.e. for /3 = a — | ^ 1 

Mla = \(l--P'^' + -,P + -P'^' + 0{P')) ■ a>a,„, . (24) 

4 V TT TT TT / 

To discuss the magnetic long range order in the whole parameter range we use instead 
of the ferro- and antiferromagnetic order parameters a parameter 

1 ^ 

^ = ]^ E (25) 
«' j'=i 

P measures the total pair correlations. For coUinear ferro- (i.e. for dominating Ji) 
or coUinear antiferromagnetic order (i.e. for dominating J2) P coincides with the 
corresponding order parameters and M^. For the LM system P is given by 

1 1 1 

8 ~ 64^ 

Because for LC and SL the order parameter P is not known in the thermodynamic 
limit we have calculated P with = 24 by exact diagonalization. The results for P 
are presented in figure 5. In accordance with the behaviour of the spin correlations 
presented in figure 4 there is a rapid change from the ferromagnetic long-range order 
for dominating Ji via a canted structure to a state with dominating antiferromagnetic 
correlations. This antiferromagnetic state is long-range ordered for the LM and SL but 
possesses no long-range order in LC because of the extremely large quantum fiuctuations 
in one dimension. There is a sharp minimum for LM and SL indicating the area with 
weak pair correlation. This minimum coincides with the maximum of the energy, i.e. 
with the point of maximal frustration. For LC the behaviour is qualitatively different, 
because in contrast to SL and LM the phase for large a is not long-range ordered. 

Of course, for SL and LC we have to take into consideration the finite size effects. 
Due to the contributions of the short range correlations P is overestimated for finite 
systems, i.e. P{N) > P{oo). For instance, the value for P larger than | at small a and 
the non-vanishing P for the LC at larger values of a is a consequence of the finite size. 
Furthermore, we can expect that for SL the minimum in P at strong frustration goes 
to zero in the thermodynamic limit, i.e. there is no long-range order for a in a region 
around 0.3. 



p = i + 

8 



a > acrtt ■ (26) 
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4.2. Thermodynamics for the Lieb-Mattis star 

The LM star is distinguished from the other systems by the existence of two additional 
integrals of motion, namely the square of the sublattice spins S^^^^^j = Y.i(zA{B) 
[Sr = Sr^ + Sjig). This fact allows to find explicitly all eigenvalues Es^r,rA,rB 

^ = + 1) - r{r + 1) - 7] + -^'Urir + 1) - r^(r^ + 1) - rB{rB + 1)] 

ziV 4 iV^ 

where we have for the relevant quantum numbers 

N 



rA(B) 



s 



0,1, 



\rA - rB\ 



for S 



J-A + r-B 



F<-A(B) ' 

for S^, 



1, 1 
2'' 2 



for Si 



The degeneracy to the quantum numbers rA(B) 



_ 2rAiB) + 1 ( 

" f + VAiB) + 1 ^ ^ 



N 
2 



(27) 

(28) 
(29) 
(30) 

(31) 



4 ^^MB), 

In addition we have to take into account the Kramers degeneracy (z-component of the 
total spin) as (2s + 1) to calculate the partition function Z 

N/4 N/4 



z 



f3J: 



exp 



rA=0 rB=0 



.2 Jo 



■fA+rB 

exp 

r = \rA-rB\ 



4J2 



Pj^[rA{rA + l) + rB{rB + 1)] 
r+1/2 



2N 



r{r + 1) 



-(3—s(s + l] 

' 2N ^ ' 



(32) 



Y (2s + 1) exp 

s=|r-l/2| 

with /3 = l/[kBT). In the thermodynamic limit the saddle-point approximation becomes 
exact, i.e. the sum in (32) is determined by its largest term. Due to symmetry we have 
= I'B- Defining the normalized sublattice polarization 



4 4 
y = -^^A = —r-B, < y < 1 



(33) 



we have for the total polarization of the medium < x < y. There exist two phases, 
which are described by the characteristics of spin correlation and the order parameters of 
the medium. One is the ferromagnetic phase (F) for dominating Ji (a < j and the other 
one is the canted (or twisted) phase (C) realized only if a > |. This latter one goes over 
smoothly in the antiferromagnetic phase for a ^ 1. The thermodynamic equations are 
as follows: 
Case 1: J2 < jJi- 

Here only the ferromagnetic phase is realized, i.e. the polarizations of the sublattice A 
and B are parallel and consequently we have r = + or x = y. The free energy per 
site is calculated as 



F = Ff = -kBTg{y) + -{hy^ - J,y) 



(34) 
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with 



2g(y) = 2 ln(2) - (1 + y) ln(l + y) - (1 - y) ln(l - y). 
The sublattice polarization y is determined via a self-consistency equation 

tanh 



y 



The spin correlations and order parameters are 

1 , 1 



\SoSi 



1 2 

4^ 



= ■ P = 



-y 



\ S,;S 



4^ 



(35) 
(36) 

(37) 
(38) 



Case 2: J2 > jJi. 

In this case we have two phases separated by a second order transition at a critical 
temperature Tq (see figure 6). Below Tq we have the canted C phase and above Tq the 
F phase described by the relations (34-38). In the C phase we have x = l/(4a) and the 
free energy reads 



F = Fc 



-ksTgiy) 



A'^^y -32X 



where the sublattice polarization y is determined by 

■1 



y 



tanh 



-.P'hy 



For the spin correlation and the order parameters we get 
1 , 11 



16a 
1 



\^i^jj\ieA,jeB 
1 o 



Mt 



32a2 
1 



-y' ; 



P 



1 



64a2 



1 2 



(39) 

(40) 

(41) 
(42) 



64a2 ' « 4^ 64a2 
The critical temperature Tq determined by 

yiT = Tg) = ]- 

4a 

is presented in figure 6. For small Tq (that means for J2 
analytic expression 

V2 8a/ 

The specific heat and the spin correlations are presented in figures 7 and 8. The second 
order transition for a > | is refiected by the kink in the spins correlations as well as in 
the peak of the specific heat. The molecular field like character of the LM system in the 
thermodynamic limit is seen by th shape of the specific heat curve (figure 7). Interesting 



(43) 

Ji/4) we have an explicit 

(44) 
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is the low temperature behaviour for strongly competing Ji and J2. Similar to the finite 
system we find order from disorder behaviour, however, this time in the correlation of 
the medium (sjSj) and not in (soSj) as for the finite system (compare figure 2 and figure 
8). In particular, for maximum frustration at a = I/a/S the correlation SiSj\i^A,]eB 
is completely suppressed for zero temperature and increases with T until the second 
order transition at Tq (figure 8). We argue that the competition between Ji and J2 is 
infiuenced by thermal fiuctuations and effectively the Ji coupling starts to overcome the 
competing J2 at finite T. For a > I/a/S the correlation function even changes its sign 
with increasing temperature. 

5. Conclusions 

In conclusion, we present the Heisenberg star with competing interactions (frustration) 
as an example for a solvable quantum spin system. The solution can be given whenever 
the solution of the surrounding medium is known. We consider as embedding media 
the linear chain, the square lattice (where only an approximative solution is given) as 
well as the Lieb-Mattis antiferromagnet which represent systems with different strengths 
of the quantum fluctuations. The effect of competing interactions manifests itself in 
a maximum of the ground state energy at maximal frustration and in a weakening of 
the magnetic correlation. In the region of maximal frustration there is a rapid change 
of correlation functions when the strength of competition is varied. This change of 
correlation is in particular dramatic when the quantum fluctuations are strong (LC). The 
competition between the exchange interactions may yield at flnite temperatures order 
from disorder phenomena, i.e. the strength of magnetic correlation can be increased by 
thermal fluctuations. 

Acknowledgments 

This work was supported by the Deutsche Forschungsgemeinschaft (project Ri 615/1-2). 
References 

[1] J.Richter and A.Voigt, J.Phys.A: Math. Gen. 27, 1139 (1994). 

[2] H.Bethe, Z. Phys. 71, 205 (1931). 

[3] R.B. Griffiths, Phys. Rev. 133, A768 (1964). 

[4] C.N. Yang and CP. Yang, Phys. Rev. 150, 321 (1966). 

[5] C.K.Majumdar and D.K.Ghosh, J. Math. Phys. 10, 1388 (1969). 

[6] I.Affieck, T.Kennedy, E.H.Lieb, and H.Tasaki, Commun. Math. Phys. 115, 477 (1988). 
[7] A.Pimpinelli, Phys. Rev. B 43, 3710 (1991). 
[8] F.D.M.Haldane, Phys. Rev. Lett. 60, 635 (1988). 
[9] B.S.Shastry, Phys. Rev. Lett. 60, 639 (1988). 



12 



[10] F.D.M.Haldane, Phys. Rev. Lett. 66, 1529 (1991). 

[11] F.D.M.Haldane, Z.N.C. Ha, J.C. Talstra, D. Bernard, and V. Pasquier, Phys. Rev. Lett. 69, 2021 
(1992). 

[12] E. Lieb and D.C. Mattis, J. Math. Phys. 3, 749 (1962). 

[13] C.Kaiser and LPeschel, J.Phys.A: Math. Gen. 22, 4257 (1989). 

[14] T.A. Kaplan, W. von der Linden and P. Horsch, Phys.Rev.B 42, 4663 (1990). 

[15] LG. Gochev and N.S. Tonchev, Phys. Rev. B 45, 480 (1992). 

[16] The Lieb-Mattis level ordering was proved rigorously for bipartite Heisenberg antiferromagnets. 
However, various numerical calculations for frustrated (non-bipartite) lattices indicate that this 
theorem holds for a large class of antiferromagnetic Heisenberg models. (see e.g. J. Richter, N.B. 
Ivanov, A.Voigt and K. Retzlaff, J. Low Temp. Phys. 99, 363 (1995) ). 

[17] J. Villain, R. Bidaux, J.-P. Caron, and R. Conte J. Phys. (Paris) 41, 1263 (1980). 

[18] C.L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 

[19] A.Chubukov, Phys. Rev. Lett. 69, 832 (1992). 

[20] J. Richter, A.Voigt, S.Kriiger, and C.Gros, Europhys.Lett. 28, 363 (1994). 

[21] A.Aharony, R.J.Birgeneau, A.Coniglio, M.A.Kastner, and H.E.Stanley, Phys. Rev. Lett. 60, 1330 
(1988); 

K. Lee and P. Schlottmann, Phys. Rev. B 42, 4426 (1990). 
[22] J. Richter , Phys. Rev. B 47, 5794 (1993). 
[23] D.D.Betts and J.Oitmaa, Phys.Rev.B 48, 10602 (1993). 
[24] J.Oitmaa, D.D.Betts and M.Aydin, Phys.Rev.B 51, 2896 (1995). 
[25] J.D.Reger and A.P.Young, Phys. Rev. B 37, 5978 (1988). 
[26] H. Neuberger and T.Ziman, Phys.Rev.B 39, 2608 (1989). 
[27] B.Bernu, C.Lluillier, and L.Pierre, Phys. Rev. Lett. 69, 2590 (1992). 
[28] E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 

[29] H. Hellmann, Einfiihrung in die Quantenchemie, Leipzig, 1937 ; R.P. Feynman, Phys. Rev. 56, 340 
(1939). 



13 



5ure captions 



^0.0 1 

Ul - 
o 

m : 
^0.1 : 



-0.2 



-0.3 : 



LC 

BCC 
LM.SL 



-0.4 



0.0 0.5 1.0 1.5 n 2.0 



Figure 1. Ground state spin correlation (sqSj) versus a = J'z/ Ji for a finite cluster 
of a central spin sq with # = 8 neighbouring spins s,- for various arrangements of 
neighbouring spins (LC - linear chain, BCC - elementary unit of a body centered 
lattice, LM - Lieb-Mattis, SL - square lattice). 




Figure 2. Spin correlations (sqSj) (a) and (sjSj) [i and j are neighbouring sites) (b) 
versus temperature T at various a = J'z/ Ji for a finite cluster of a central spin sq with 
# = 8 neighbouring spins s,- which form a linear chain (LC). 
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Figure 3. Ground state energy versus frustration parameter a for N ^ co. 
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Figure 4. Ground-state spin correlation within the medium (sjSj) and between the 
central spin and the embedding medium (sqSj) versus a in the thermodynamic limit 

# CXD. 
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Figure 5. Ground-state total correlation parameter P (cf. equation (25)) of the 
medium versus a. 




Figure 6. Transition temperature Tq versus a for the LM system and illustration of 
the principal arrangement of the spins in the ferromagnetic F phase and the canted C 
phase. The lower arrow represents the central spin and the two upper arrows represent 
the sublattice spins of the LM system. 
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Figure 7. Specific heat of the LM star versus temperature T for several a. 



